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ABSTRACT 

We present precise measurements of the X-ray gas mass fraction for a sample of lumi- 
nous, relatively relaxed clusters of galaxies observed with the Chandra Observatory, 
for which independent confirmation of the mass results is available from gravitational 
lensing studies. Parameterizing the total (luminous plus dark matter) mass profiles 
using the model of Navarro, Frenk & White (1997), we show that the X-ray gas mass 
fractions in the clusters asymptote towards an approximately constant value at a ra- 
dius r2500j where the mean interior density is 2500 times the critical density of the 
Universe at the redshifts of the clusters. Combining the Chandra results on the X-ray 
gas mass fraction and its apparent redshift dependence with recent measurements of 
the mean baryonic matter density in the Universe and the Hubble Constant deter- 
mined from the Hubble Key Project, we obtain a tight constraint on the mean total 
matter density of the Universe, f2 m = 0.30+q"q3, and measure a positive cosmological 

constant, Q\ = 0.951q^. Our results are in good agreement with recent, independent 
findings based on analyses of anisotropies in the cosmic microwave background radia- 
tion, the properties of distant supernovae, and the large-scale distribution of galaxies. 
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1 INTRODUCTION 

The matter content of rich clusters of galaxies is thought to 
provide a fair sample of the matter content of the Universe as 
a whole (White et al. 1993). The observed ratio of the bary- 
onic to total mass in clusters should therefore closely match 
the ratio of the cosmological parameters fib/f^m, where fib 
and fi m are the mean baryon and total mass densities of the 
Universe, in units of the critical density. The combination of 
robust measurements of the baryonic mass fraction in clus- 
ters with accurate determinations of from cosmic nucle- 
osynthesis calculations (constrained by the observed abun- 
dances of light elements at high redshifts) can therefore be 
used to determine S7 m . 

This method for measuring Q m , which is particularly 
simple in terms of its underlying assumptions, was first high- 
lighted by White & Frenk (1991) and subsequently employed 
by a number of groups (e.g. Fabian 1991, White et al. 1993, 
David, Jones & Forman 1995; White & Fabian 1995; Evrard 
1997; Fukugita, Hogan & Peebles 1998; Ettori & Fabian 
1999; Bahcall et al. 2000). In general, these studies have 
found fi m < 1 at high significance, with preferred values 
lying in the range fi m ~ (0.1 - 0.3)^~ 5 . 

Sasaki (1996) and Pen (1997) described how measure- 
ments of the mean baryonic mass fraction in clusters as a 
function of redshift can, in principle, be used to place more 



detailed constraints on cosmological parameters, since the 
observed baryonic mass fraction values are sensitive to the 
angular diameter distances to the clusters assumed in the 
analyses. Until now, however, systematic uncertainties in the 
observed quantities have seriously complicated the applica- 
tion of such methods. 

The baryonic mass content of rich clusters of galaxies is 
dominated by the X-ray emitting intracluster gas, the mass 
of which exceeds the mass of the optically luminous material 
by a factor ~ 6 (e.g. White et al. 1993; David et al. 1995; 
Fukugita, Hogan & Peebles 1998) . Since the X-ray emissivity 
of the X-ray gas is proportional to the square of its density, 
the gas mass profile can be precisely determined from the 
X-ray data. With the advent of accurate measurements of 
Qh (e.g. O'Meara et al. 2001 and references therein) and 
a precise determination of the Hubble Constant (Freedman 
et al. 2001), the dominant uncertainty in determining Q m 
from the baryonic mass fraction in clusters has lain in the 
measurements of the total (luminous plus dark) matter dis- 
tributions in the individual clusters. 

In this letter we report precise measurements of the X- 
ray gas mass fraction for a sample of luminous, relatively 
relaxed clusters spanning the redshift range 0.1 < z < 0.5, 
for which precise, consistent mass models have recently 
been determined from Chandra X-ray data and independent 
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Table 1. Summary of the Chandra observations. 





z 


Date 


Exposure (ks) 


PKS0745-191 


0.103 


2001 


Jun 


16 


17.9 


Abell 2390 


0.230 


1999 


Nov 


07 


9.1 


Abell 1835 


0.252 


1999 


Dec 


12 


19.6 


MS2137-2353 


0.313 


1999 


Nov 


18 


20.6 


RXJ1347-1145(1) 


0.451 


2000 


Mar 


05 


8.9 


RXJ1347-1145(2) 


0.451 


2000 


Apr 


29 


10.0 


3C295 


0.461 


1999 


Aug 


30 


17.0 



gravitational lensing constraints (Allen et al. 2001a, 2002; 
Schmidt, Allen & Fabian 2001). The agreement between the 
mass measurements determined from the two independent 
methods firmly limits the systematic uncertainties in the 
baryonic mass fraction measurements to < 10 per cent, an 
accuracy comparable to the current flb results. With the 
reduced systematic uncertainties, we show that a method 
similar to those proposed by Sasaki (1996) and Pen (1997) 
can be successfully applied to the data, resulting in a tight 
constraint on n m and an interesting constraint on Qa- We 
show that the results obtained are in good agreement with 
those from recent studies of anisotropies in the cosmic mi- 
crowave background radiation, the large-scale distribution 
of galaxies, and the properties of distant supernovae (e.g. 
Jaffe et al. 2001; Efstathiou et al. 2001). 

Results on the X-ray gas mass fractions are 
quoted for two default cosmologies: SCDM with h = 
Ho/100 km s" 1 Mpc _1 = 0.5, Sl m = 1 and fi A = 0, and 
ACDM with h = 0.7, fl m = 0.3 and f2 A = 0.7. 



2 OBSERVATIONS AND DATA ANALYSIS 

The Chandra observations were carried out using the back- 
illuminated S3 detector on the Advanced CCD Imaging 
Spectrometer (ACIS) between 1999 August 30 and 2001 
June 16. For our analysis we have used the the level-2 event 
lists provided by the standard Chandra pipeline processing. 
These lists were cleaned for periods of background flaring us- 
ing the CIAO software package resulting in the net exposure 
times summarized in Table |l| 

The Chandra data have been analysed using the meth- 
ods described by Allen et al. (2001a, 2002) and Schmidt 
et al. (2001). In brief, concentric annular spectra were ex- 
tracted from the cleaned event lists, centred on the peaks 
of the X-ray emission from the clusters.^ The spectra were 
analysed using XSPEC (version 11.0: Arnaud 1996), the 
MEKAL plasma emission code (Kaastra & Mewe 1993; in- 
corporating the Fe-L calculations of Liedhal, Osterheld & 
Goldstein 1995), and the photoelectric absorption models of 
Balucinska-Church & McCammon (1992). Only data in the 
0.5 — 7.0 keV energy range were used. The spectra for all 
annuli were modelled simultaneously, in order to determine 



For RXJ1347-1145, the data from the southeast quadrant of 
the cluster were excluded due to ongoing merger activity in that 
region; Allen et al. (2002). 



the deprojected X-ray gas temperature profiles under the 
assumption of spherical symmetry. 

For the mass modelling, azimuthally-averaged surface 
brightness profiles were constructed from background sub- 
tracted, flat-fielded images with a 0.984 x 0.984 arcsec 2 pixel 
scale (2x2 raw detector pixels) . When combined with the de- 
projected spectral temperature profiles, the surface bright- 
ness profiles can be used to determine the X-ray gas mass 
profiles (to high precision) and total mass profiles in the clus- 
ters.]^ For this analysis, we have used an enhanced version 
of the image deprojection code described by White, Jones 
& Forman (1997) with distances calculated using the code 
of Kayser, Helbig & Schramm (1997). 

We have parameterized the cluster mass (luminous plus 
dark matter) profiles using a Navarro, Frenk & White (1997; 
hereafter NFW) model with 



p(r) 



p c (z)S c 



(r/r s ) (1 + r/r s y 



(1) 



where p(r) is the mass density, p c (z) = 3H(z) 2 /8ivG is the 
critical density for closure at redshift z, r s is the scale ra- 
dius, c is the concentration parameter (with c = r2oo/r s ) 
and 5 C = 200c 3 /3 [ln(l + c) - c/(l + c)]. The normalizations 
of the mass profiles may also be expressed in terms of an 
equivalent velocity dispersion, a — \ // 50r B cH(z) (with r B in 
units of Mpc and H(z) in kms Mpc" 1 ). 

In determining the results on the X-ray gas mass frac- 
tion, / gas , we have adopted a canonical radius r2soo, within 
which the mean mass density is 2500 times the critical den- 
sity of the Universe at the redshift of the cluster. (The T2500 
values are determined directly from the Chandra data, with 
confidence limits calculated from the \ 2 grids.) The r2soo 
values are well-matched to the outermost radii at which 
reliable temperature measurements can be made from the 
Chandra data. Note that the data for PKS0745-191 do not 
quite reach to r25oo and for this cluster we quote / gas at 
the outermost radius at which reliable measurements can 
be made ~ 0.8r25oo- Although independent confirmation of 
the X-ray mass results for 3C295 is not available, we in- 
clude this cluster in our simple since in most other ways 
it appears similar to the other objects in the sample. The 
?"2500 values for the clusters are listed in Table ^[ The best-fit 
NFW model parameters and 68 per cent confidence limits 
are summarized by Allen et al. (2001b). 



T The observed surface brightness profile and a particular pa- 
rameterized mass model are together used to predict the temper- 
ature profile of the X-ray gas. (We use the median temperature 
profile determined from 100 Monte-Carlo simulations. The outer- 
most pressure is fixed using an iterative technique which ensures 
a smooth pressure gradient in these regions.) The predicted tem- 
perature profile is rebinned to the same binning as the spectral 
results and the \ 2 difference between the observed and predicted, 
deprojected temperature profiles is calculated. The parameters for 
the mass model are stepped through a regular grid of values in 
the r B -cr plane (see text) to determine the best-fit values and 68 
per cent confidence limits. (The best-fit models generally provide 
good descriptions of the data). The gas mass profile is determined 
to high precision at each grid point directly from the observed sur- 
face brightness profile and model temperature profile. Spherical 
symmetry and hydrostatic equilibrium are assumed throughout. 
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2500 



Figure 1. The observed X-ray gas mass fraction profiles with 
the radial axis scaled in units of r2500- Symbols are as follows: 
PKS0745-191 (light circles), Abell 2390 (light triangles), Abell 
1835 (dark triangles), MS2137-2353 (light squares), RXJ1347- 
1145 (dark circles), 3C295 (dark squares). The default ACDM 
cosmology is assumed. Note that / gas (r) is an integrated quan- 
tity and so the error bars on neighbouring points in a profile are 
correlated. 



3 RESULTS 

3.1 The X-ray gas mass fraction measurements 

Fig. [I] shows the observed f ga s(r) profiles for the six clus- 
ters assuming the standard ACDM cosmology. We see that 
although some variation is present from cluster to cluster, 
the profiles tend towards a similar value at r25oo- Table B 
lists the results on the X-ray gas mass fractions measured 
at r 2 500 for both the SCDM and ACDM cosmologies. Tak- 
ing the weighted mean of the / gas results for all six clusters 
studied, we obtain / gas = 0.160 ± 0.007 for SCDM (h=0.5) 
and / gas = 0.113 ± 0.005 for ACDM (h=0.7). 

In calculating the total baryonic mass in the clusters, we 
assume that the optically luminous baryonic mass in galaxies 
is 0.19/i 0,5 times the X-ray gas mass (White et al. 1993; 
Fukugita, Hogan & Peebles 1998) . Other sources of baryonic 
matter are expected to make only very small contributions 
to the total mass and are ignored. 

Given the baryonic masses, and assuming that the re- 
gions of the clusters within r2soo provide a fair sample of the 
matter content of the Universe, we can write 



a, 



fi h 



/gas(l + 0.19/1 "' 



(2) 



For Q, h h 2 = 0.0205 ± 0.0018 (O'Meara et al. 2001) and us- 
ing the ACDM (h = 0.7) / gas values, we obtain the (self- 
consistent) result fi m = 0.319 ± 0.032. Using the SCDM 
(h = 0.5) / gas values, we obtain fi m = 0.452 ± 0.044. 



3.2 Cosmological constraints from the f gas (z) data 

In addition to the simple calculation of fi m based on the 
weighted-mean / gas values, described above, the data for 
the present sample can be used to obtain more rigorous 



constraints on cosmological parameters from the apparent 
variation of / gas with redshift. 

Fig. ^ shows the / gas values as a function of redshift for 
the SCDM and ACDM cosmologies. We see that whereas 
the results for the ACDM cosmology are consistent with a 
constant / gas value, the results for SCDM indicate an ap- 
parent drop in / gas as the redshift increases. The differences 
in the f gas (z) behaviour for the SCDM and ACDM cos- 
mologies reflect the dependence of the / gas (z) measurements 
on the assumed angular diameter distances to the clusters 
(/ gas oc Da' 5 )- Under the assumption that the / gas values 
should be invariant with redshift, as would be expected if 
rich, relaxed clusters provide a fair sample of the matter 
content of the Universe, we can see from inspection of Fig. ^ 
that the data for the present sample favour the ACDM over 
the SCDM cosmology. 

In order to quantify more precisely the degree to which 
our data can constrain the relevant cosmological param- 
eters, we have fitted the data in Fig. |i](a) with a model 
which accounts for the expected apparent variation in the 
fgas(z) values, which are measured assuming an SCDM cos- 
mology, for different underlying cosmologies. The 'true' cos- 
mology should be the cosmology that provides the best fit 
to the measurements. (We work with the SCDM data. Note 
that the / gas (r) profiles exhibit only small variations around 
f 2500 , and so the effects of changes in r25oo as the cosmology 
is varied can be ignored.) 

The model function fitted to the data is 



/•mod, s 
J gas \^ I 



fib 



i + o.wVK) fi„ 



05 



=i,n A =o 



D 



(3) 



which depends on fi m , fiA, fib and h. The ratio (Zi/0.5) 1 ' 
accounts for the change in the Hubble Constant between the 
considered model and default SCDM cosmology, and the ra- 
tio of the angular diameter distances accounts for deviations 
in the geometry of the Universe from the Einstein-de Sitter 
case. We constrain Q b h 2 = 0.0205 ± 0.0018 (O'Meara et al. 
2001) and h = 0.72 ± 0.08, the final result from the Hub- 
ble Key Project reported by Freedman et al. (2001). The x 2 
difference between the model and SCDM data is then 



E 



[/gas ( z i) fg. 



all clusters 



' f 



+ 



fib/i 2 - 0.0205 \ 2 (h- 0.72 



0.0018 



0.08 



(4) 



where / gas , ; and ff/ gaB ; are the best-fit values and symmetric 
root-mean-square errors for the SCDM data from Table H, 
and Zi are the redshifts of the clusters. We have examined a 
grid of cosmologies covering the plane 0.0 < fi m < 1.0 and 
0.0 < fiA < 1-5. The joint 1, 2 and 3 a confidence contours 
on fi m and fiA (corresponding to A% 2 values of 2.30, 6.17 
and 11.8, respectively) determined from the fits are shown 
in Fig. |. 

The best-fit cosmological parameters and marginalized 
Iff error bars are fi m = 0.301^ and fi A = 0.95t ;^, 
with Xmin = 1-7 for 4 degrees of freedom, indicating that 
the model provides an acceptable description of the data. 
The best-fit cosmological parameters are similar to those 
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Table 2. The observed X-ray gas mass fractions (and 
and ACDM cosmologies. 



per cent confidence limits) measured at r2500 (in Mpc) for the default SCDM 



SCDM 



ACDM 



''■jr,(io 



''2r,()(l 



PKS0745-191 
Abell 2390 
Abell 1835 
MS2137-2353 
RXJ1347-1145 
3C295 



0.69 4 



U - S5 -0.05 
3+0.14 
-0.09 
n 79 +0.05 
U -' Z -0.03 
n 49+0-03 
U -* M -0.01 

u -' -0.08 
42+ - 03 



174+0-013 
u - 1,y: -0.012 

209+ 060 
u.zu»_ 046 

164+0-016 

159+ 009 
q 142 +0.034 
u - 1 ^-0.027 
n 1OO+0 020 
U - ii8 -0.016 



n fis+0 03 

,+0.15 
-0.09 
q+0.06 
-0.02 
q +0.02 
-0.03 
n 7 o+0.08 
u - lo -0.09 
41+0- 04 



0.64] 
0.66^ 
0.46"* 



n i -, 9 +0.008 
u - xl -0.009 
n 1,0+0.047 
u.loo_ 033 

H4+0.006 
u ' llyi -0.013 
-,-,-+0.015 
v - 11 1 -0.009 

n ins+0 03i 
u.iuo_ 016 




0.5 




0.5 



Figure 2. The apparent variation of the observed X-ray gas mass fraction (with root-mean-square la errors) as a function of redshift for 
the default (a: left panel) SCDM (h=0.5) and (b: right panel) ACDM (h=0.7) cosmologies. The dashed curves show the results of fitting 
a constant value to the data in each case. The solid line in (a) shows the predicted curve for the best-fit cosmology with fi m = 0.30 and 
f2 A = 0.95 (see Section 3.2). 



assumed for the default ACDM cosmology in Fig. gb, which 
is expected given the approximately constant nature of the 
/gas(z) values shown in that Figure. 



4 DISCUSSION 

The result on the mean matter density of the Universe, 
Qm = 0.301q 03> determined from the Chandra results on the 
X-ray gas mass fraction for the present sample of relaxed, 
lensing clusters, represents one of the tightest constraints 
on this cosmological parameter to date. The variation of the 
gas mass fraction with redshift also yields the measurement 
of a positive cosmological constant with Qa = 0.951q'^, 
in good agreement with previous results based on studies 
of the properties of distant supernovae (Riess et al. 1998; 
Perlmutter et al. 1999) 

In Fig. ^ we show a comparison of the joint con- 
straints on f2 m and Qa determined from the Chandra / gas (z) 
data, with the results of Jaffe et al. (2001) from studies 
of cosmic microwave background (CMB) anisotropies (in- 
corporating the COBE Differential Microwave Radiometer, 



BOOMERANG-98 and MAXIMA- 1 data of Bennett et al. 
1996, de Bernardis et al. 2000 and Hanany et al. 2000, re- 
spectively)^] and the properties of distant supernovae (in- 
corporating the data of Riess et al. 1998 and Perlmutter et 
al. 1999). The agreement between the results obtained from 
the independent methods is striking: all three data sets are 
consistent, at the la confidence level, with a cosmological 
model with f2 m = 0.3 and Qa = 0.7 — 0.8. These results are 
also consistent with the findings of Efstathiou et al. (2001) 
from a combined analysis of the 2dF Galaxy Redshift Survey 
and CMB data. 

An important aspect of the present work is that, in 
addition to the exquisite data quality provided by Chan- 
dra, the clusters studied are all regular, relatively relaxed 
systems for which independent confirmation of the mass re- 



t We note that the results on f2 m and £7a from the CMB data 
reported by Jaffe et al. (2001) are consistent with, though less 
constraining than, the later analyses of de Bernardis et al. (2002) 
and Stompor et al. (2001) using the full BOOMERANG and 
MAXIMA-1 data sets. 
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suits is available from gravitational lensing studies. The sys- 
tematic uncertainties in the / ga s measurements are therefore 
greatly reduced with respect to previous X-ray studies. For 
both Abell 2390 and RXJ1347-1145, the X-ray and weak 
lensing mass profiles are consistent within their 68 per cent 
confidence limits. For Abell 1835, 2390, MS2137-2353 and 
PKS0745-191, the observed strong lensing configurations (on 
scales r ~ 20 — 80 h~ kpc) can be explained by mass mod- 
els within the 68 per cent Chandra confidence contours, al- 
though redshift measurements for the arcs (which are re- 
quired to define the lensing masses precisely) are not avail- 
able in all cases ]§] The presence of significant non-thermal 
pressure support (e.g. arising from turbulent and/or bulk 
motions and/or magnetic fields) on scales ~ r25oo can there- 
fore be excluded, and the residual systematic uncertainties 
in the / gas values are small (< 10 per cent i.e. smaller, typ- 
ically, than the statistical uncertainties. We note that the 
effects of departures from spherical symmetry on the / gas 
results are expected to be < a few per cent e.g. Buote & 
Canizares 1996). 

The observed f sas (r) profiles are essentially flat around 
^2500, which supports the assumption that the measured / gas 
values represent a fair sample of the matter content of the 



§ For RXJ1347-1145, a two-component mass model, consistent 
with the complex X-ray structure observed in the southeast quad- 
rant, is required to explain the strong lensing data. 



Universe. If, however, the values were to rise by a further 
~ 10 per cent beyond r2soo, the result on Q m would drop 
by a corresponding amount. The / gas values measured at 
T2500 are not sensitive to the choice of using an NFW model 
to parameterize the total mass distributions in the clusters. 
Repeating the analysis presented here using either a non- 
singular isothermal sphere or a Moore et al. (1998) model to 
parameterize the total mass distributions leads to results on 
the weighted-mean / gas values in good agreement with those 
quoted in Section 3.1. We note, however, that the f g as(r) 
profiles determined using the different mass models exhibit 
some systematic variation, particularly at small radii (r < 
0.1r25oo), and when extrapolated to large radii (r > r2soo), 
as can be expected given the different asymptotic slopes. 
In a future paper, we will examine in detail the ability of 
different parameterized models to describe the Chandra data 
for relaxed clusters. 

The constraints on f2 m and Qa should improve as fur- 
ther Chandra, XMM-Newton and high-quality gravitational 
lensing data become available for more regular, relaxed clus- 
ters, especially at high redshifts (although relaxed systems, 
like those studied here, are expected to be very rare at high 
redshifts). This work can also be extended to include less 
relaxed clusters, or clusters which appear relaxed at X-ray 
wavelengths but for which independent confirmation of the 
mass results from gravitational lensing studies is not yet 
available, although this will require careful consideration of 
the additional systematic uncertainties involved. 
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